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We present a new type of calorimeter in which we couple an unknown heat capacity 
with the aid of Peltier elements to an electrical circuit. The use of an electrical in- 
ductance and an amplifier in the circuit allows us to achieve autonomous oscillations, 
and the measurement of the corresponding resonance frequency makes it possible to 
^ I accurately measure the heat capacity with an intrinsic statistical uncertainty that 

'^ ■ decreases as ~ tm with measuring time tm, as opposed to a corresponding uncer- 

tainty ~ tm in the conventional alternating current (a.c.) method to measure heat 
capacities. We have built a demonstration experiment to show the feasibility of the 
*7i • new technique, and we have tested it on a gadolinium sample at its transition to the 

ferromagnetic state. 
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I. INTRODUCTION 

Calorimetry is an old field in experimental physics, and numerous experimental techniques 
have been developed for the accurate measurement of heat capacities^. The frequently used 
a.c. technique invented by Corbino^ and modified by Sullivan and Seidell, for example, mea- 
sures a heat capacity by monitoring the temperature amplitude of a sample that is subject 
to an oscillating heating power , or, although less common, by measuring an oscillating heat 
current produced by a Peltier element^. Such a.c. techniques have numerous advantages 
over other methods, e.g., heat-pulse, differential-thermal-analysis^ or relaxation techniques^. 
A certain robustness to noise and the possibility to use very small temperature amplitudes 
stems from the fact that the signal of interest can be separated from the instrumental noise 
by using a frequency selective filter, such as a lock-in amplifier. To increase the total ac- 
curacy even further, the duration of a measurement tm can be increased at will, and the 
resulting statistical uncertainty decreases as tm . 

If a calorimeter could be built where the heat capacity is related to the frequency instead of 
the amplitude of an oscillation one would retain most of the advantages of the conventional 
a.c. techniques, but the statistical uncertainty would decrease much more rapidly^, i.e., as 
tm • Therefore, the accuracy of such a method would be mainly limited by the stability of 
the oscillator components, rather than by statistical constraints. 

In the following we describe the realization of an oscillating thermo-electrical circuit using 
Peltier elements which makes an electrical inductor to effectively act as a "thermal induc- 
tor". The heat capacity can be calculated from the measured resonance frequency of the 
circuit, and its statistical uncertainty drops according to tm as expected. 



II. OPERATING PRINCIPLE OF THE THERMO-ELECTRICAL 
OSCILLATOR 

It is known that the flow of heat Qt (with a corresponding heat current 1% = Qt) between 
entities with a certain heat capacity Ct through links with a finite thermal conductivity can 
be mapped onto the flow of charge Qe (with the electrical current Jc = Qc) in electrical 
circuits that is composed of capacitors C^ and resistors, respectively, where the tempera- 
ture T in the thermal case corresponds to the electrostatic potential V in the equivalent 




FIG. 1. An electrical inductor Lc is connected to a Peltier element that separates a heat capacity 
Ct from a thermal bath T. Changes in the electrical current Iq induce a temperature difference 
AT = Lc/c/a (where a denotes the Seebeck coefficient) across the Peltier element that counteracts 
the variation in the flow of heat It- 



electrical circuit. The differential equations describing the time evolution of Qt (t) or Qc (t), 
respectively, are of first order in time t. However, a simple thermal analogue of an elec- 
trical inductor involving a higher derivative Qt = -^t does not exist, although the inclusion 
of terms including It is justified in certain rare cases^. Therefore it is impossible to realize 
autonomously oscillating thermal circuits by using only standard thermal elements. 

In our experiment we connect an electrical inductor Lc to a Peltier element as shown in 
Fig. [U If the heat current It through the Peltier element to the thermal bath varies with 
time, an associated electrical current /e in the circuit changes according to IIJo = It , where 
n is the Peltier coefficient of the Peltier element. It induces a voltage drop L^Ic across the 
inductor that opposes the changes in the heat current It through the Peltier element via 
the Seebeck effect, i.e., a resulting additional temperature difference AT = Lch/oi (where a 
denotes the Seebeck coefficient) across the Peltier element is imposed that counteracts the 
variation in the flow of heat. It ■ We can therefore define an effective thermal inductance 
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However, as a consequence of the significant losses produced by the non-zero thermal con- 
ductivity and electrical resistance of the Peltier element, the idealized element shown in 



Fig. [Hallows in general only for an overdamped variation of the heat current. This obsta- 
cle can be avoided by including an amplifier with a suitable gain in the electrical circuit, 
and autonomous oscillations can eventually be achieved in this way. In the ideal case of a 
stationary oscillation with constant amplitude where all losses are exactly compensated by 
the amplifying element, the electrical circuit behaves like a simple LC circuit with thermal 
inductance L^ and a heat capacity Ct, or with an electrical inductance L^ and an equivalent 
electrical capacitance 

"" V a AT a/n/edt aU ^ ' 

and with resonance frequency u"^ = 1/LtCt = 1/ L^Cc- Therefore, the heat capacity can be 
derived from uj according to 



C. = ^-,. (3) 
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notably without any need to quantify the thermal and electrical losses in the Peltier element. 
To the best of our knowledge, the correspondences between an L^ and a thermal analogue 
Lt, or Ct and a C^ as expressed in Eqs. ([T]) and ([2]), respectively, are novel and have also 
never been exploited in any calorimetric technique. 

III. EXPERIMENTAL SETUP 

We have chosen the circuit shown in Fig. [2] containing two Peltier elements (MPC-D701, 
Micropelt Inc.) based on Bi2Te3, with lateral dimensions of ~ 3.5 mm and a thickness 
of ~ 1 mm. While the left element acts, together with the sample to be measured, as the 
equivalent electrical capacitance as described above, the right element is used to compensate 
for thermal and electrical losses and to maintain a constant amplitude of the temperature 
oscillation. For this purpose, the voltage across the right element is amplified and then 
fed onto the electrical contacts of the left element. To achieve a stable oscillation of the 
desired amplitude, the gain of the amplifier can be controlled using a voltage-controlled gain 
amplifier (AD630, Analog Devices). An additional clamp with Schottky diodes (BAT43) 
and two potentiometers forms a nonlinear resistor that reduces the gain at large oscillation 
amplitudes and therefore forces the circuit to self-stabilize. 

For a design heat capacity Ct ~ 20 mJ/K to be measured near room temperature (which 
corresponds to 2.87 ■ 10~^mol of a Dulong-Petit solid), and using a ~ 29mV/K and 11 ^ 
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FIG. 2. Circuit diagram showing the Peltier element system, the Gyrator to mimic a large electrical 
inductance, and the amplifier to compensate for thermal and electrical losses (see text). The voltage 
Vg controls the gain of the amplifier. The two Peltier elements have the same polarity, and the 
sample is attached to their "top side". The measured signal is tapped at Vout- 




FIG. 3. Thermal oscillations of the LC-circuit. The solid line is a least-squares fit to a sinusoidal 
variation of AT in time around an offset value that is caused by lost heat dissipated in the Peltier 
elements. 



aT ^ 8.4 W/A for the used Peltier elements, the equivalent electrical capacitance becomes 
Co ~ 82 mF. Depending on the heat capacities and the thermal conductivities involved 
in thermal part of the system, the frequency of the temperature oscillation must be small 
enough to allow for a quasi- stationary thermal equilibrium. If we choose / = Ci;/27r ^ 0.18 Hz 
we need an inductance Lc ~ lOH which would require impractically large passive elements. 
Therefore we replaced the inductance by a gyrator, the circuit diagram of which is also shown 
in Fig. 121 The effective induction of the gyrator is given by^ L^ = {Rp + RGi)RG2CgyT, where 
the internal resistance Rp ~ 6.5 f2 of the Peltier element has to be included here. With 
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FIG. 4. a) Measured frequency of the oscillator and b) heat capacity of a Gd sample, as functions 
of temperature. Crosses represent the data that we obtained by the present thermal LC-method, 
while circles denote the values from a measurements on the same sample using a commercial heat 
capacity system (see text). 

Rqi = Ifi , Rg2 = 200 kfi and Cgyr = 6.8 /xF we indeed obtain Le ~ 10 H. 



IV. TEST MEASUREMENTS 

A test sample (45.2 nig or 2.87- 10~^ niol of gadolinium with 99.99% purity that had been 
previously annealed at T = 1123 K in vacuum to reduce lattice strains^) was sandwiched 
between the two Peltier elements and placed into a commercial PPMS system (Physical 
Property Measurement System, Quantum Design) for control of the base temperature. The 
gain of the amplifier was chosen in a way to achieve temperature oscillations with an ampli- 
tude of the order of 0.2 K. For each data point, the oscillating voltage was read out during 
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1200 s with a 10 Hz-triggered Keithley 2001 Multimeter in equal time intervals of 0.1 s 
(see Fig. |3]), and the resulting data were then least-squares fitted to obtain the oscillation 
frequency (see Fig. H^). To compensate for the heat capacity of the empty system (^ 15% 
of the total heat capacity), corresponding measurements have been done without the sample 
installed. The resulting values of the heat capacity are shown in Fig. |3]d, together with heat- 
capacity data obtained on the same sample using the commercial heat-capacity option of the 
PPMS system and with temperature steps of ~ 0.5 K. While the heat-capacity data using 
our thermal LC- method are neither smoothed nor further averaged, the PPMS data were 
averaged over 6 individual measurements. The data-acquisition time for each of the LC- 
method data points was 1200 s, while collecting an averaged PPMS data point took about 
one hour of time, and the temperature amplitudes for both techniques were comparable (0.3 
K - 0.5 K). The transition of the Gd sample to the ferromagnetic state takes place at the 
Curie temperature Tq ^ 290K, and the absolute values for the heat capacity agree very well 
with the PPMS data (see Fig. Eb). 

V. ACCURACY CONSIDERATIONS 

The accuracy of any heat-capacity measurement is limited by systematic and statistical 
uncertainties. Sources of systematic errors (such as finite instrumental resolution, limited 
accuracy by a large background signal, errors in the temperature calibration etc.) are 
common to all techniques for measuring a heat capacity, and we shall not further discuss 
these shortcomings here. The main advantage of the present thermal LC-circuit technique 
is that it significantly reduces the statistical uncertainty with increasing measuring time. 

The statistics behind the problem of estimating the amplitude A, the frequency u and 
the phase ip of an oscillation A sin(cjt -|- ip) out of a set of A^ samples of finite length T has 
been analyzed in detail in Ref. [g. Assuming a white noise with variance cr^, the resulting 
relative standard deviation a a/ A in the amplitude A from a measurement of A^ oscillation 
cycles with length T = 27r/cj is 

^ > ^ ■ (4) 
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In a conventional a.c. technique to measure heat capacities, the relative standard devia- 
tion in the heat capacity Ct from the mean value is proportional to the respective standard 
deviation in the measured amplitude A of a temperature oscillation^, and therefore the sta- 
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FIG. 5. a) Amplitudes and b) frequencies for a selected data point together with the respective 
standard deviations from the Levenberg-Marquardt fit, as functions of measurement time tm- c) 
The respective relative standard deviations on a log-log scale to illustrate the power-law behavior 
discussed in the text Solid lines are fits according to t^^'"^ for (ta/A and to t~'^' ^ for the Ctj/a; data, 
respectively, d) The overlapping Allan deviation ay, plotted against the observation time r. The 
solid line illustrates the approximate r~^'^ dependence of ay (see text). 
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tistical uncertainty in C^ decreases with tm where t^ = NT is the total time to take such 
a measurement (in all other conventional methods, the statistical uncertainty also decreases 
as A^~^/^, where A^ is the number of measurements to calculate the average of a data point). 
The corresponding relative standard deviation a^/u in the frequency u , however, varies 
with A^ according to 
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w ■nA^N{N'^ - 1) ' 
i.e., it decreases with N~^^'^ (or tm ) for large N. As the heat capacity in our technique is 



proportional to u;~^ (see Eq. E]), the relative standard deviation oc^jCx, ~ 2 ■ o^juj also de- 
creases with N in the same fashion. In Fig. |5]we show the amplitudes and the frequencies for 
a selected data point, together with the respective relative standard deviations, as functions 
of the measurement time t^ which is proportional to the number of oscillation cycles A^. 
The standard deviations a a and a^^ have been taken from a standard Levenberg-Marquardt 
routine used for least-squares fitting the experimental data. The (JaI-^ data decrease rela- 
tively slowly with 'N~^l'^ , while the o^joj data (and along with them ac^/Ci ) drop much 
faster with N''^^"^ as expected. Therefore, very high accuracies for Ct can, in principle, be 
achieved in a comparably very short measuring time. For example, to reduce the statisti- 
cal uncertainty in Ct by a factor of 10, conventional methods require a 100 times longer 
duration of a measurement, while the thermal LC method needs only a ~ 5 times longer 
data acquisition time. As both cr^/^ and (Ji^/u are proportional to the same quantity a /A, 
there is no hidden drawback in the equation [S] that would put this tremendous statistical 
advantage into perspective. 

However, in order to tap the full potential of this technique it is essential that the os- 
cillation frequency u is sufficiently stable during the time of the measurement so that the 
potentially minuscule statistical uncertainty is not impaired by extrinsic fluctuations, e.g., 
by an insufficient stabilization of the temperature of the sample or by an unwanted variabil- 
ity in the electronic circuit. A standard measure for the stability of an oscillator is provided 
by the Allan deviation (Jy—. Assuming white noise in the oscillator frequency w, for exam- 
ple, ay is expected to drop as r~^/^, where r denotes the time over which the frequency i 
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observed. At large observation times ay can saturate or even increase again as soon as other 
sources of noise (Flicker noise, random-walk noise) become dominantiii. In Fig. [5li we have 
plotted the overlapping Allan deviation as a function of observation time r on a log-log scale. 
While ay initially drops approximately with r~^/^, is saturates in our experiment around 
ay ^ A ■ 10"^ for r ~ 1000 s an increases for larger vales of r. In this sense we interpret the 
minimum in ay as an estimate for the achievable uncertainty in u for a given set of experi- 
mental data, while the a^^/uj values from Eq. [5] and shown in Fig. [St represent the statistical 
uncertainty assuming an infinitely stable frequency with white noise only in the temperature 
amplitude. In our demonstration experiment, no special care has been taken to minimize 
extrinsic sources for variations in uj. Such sources of error are the decisive factors that are 
hmiting the maximum achievable accuracy in the present thermal LC technique but they 



are to a large extent controllable, while the intrinsic statistical uncertainty should play a 
significantly less important role than in conventional techniques to measure heat capacities. 

VI. CONCLUSIONS 

We have described a method where a heat capacity is coupled to an electrical circuit 
with the aid of Peltier elements, which makes an electrical inductor to effectively act as a 
"thermal inductor". The use of an amplifier in the circuit allows us to achieve autonomous 
oscillations, and the measurement of the corresponding resonance frequency makes it possible 
to calculate the heat capacity. While the statistical uncertainty in conventional techniques 
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to measure heat capacities scales according to tm with measuring time tm, the respective 
uncertainty in the present thermal LC technique decreases significantly faster, i.e., with 
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tm , and is ultimately limited by the stability of the thermal and electronic components 
used in the experiment . 
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